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We experimentally determine the quantum discord present in two-mode squeezed vacuum gener- 
ated by a four- wave mixing process in hot rubidium vapor. The frequency spectra of the discord, 
as well as the quantum and classical mutual information are also measured. In addition, the effects 
of symmetric attenuation introduced into both modes of the squeezed vacuum on the discord, the 
quantum mutual information and the classical correlations are examined experimentally. Finally, 
we show that due to the multi-spatial-mode nature of the four- wave mixing process, the quantum 
discord may exhibit sub- or superadditivity depending on which spatial channels are selected. 
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Quantum communication and computation schemes 
have been thought to rely on entanglement to offer im- 
provements over classical systems. Entanglement as a 
necessary resource can be limiting for experimental im- 
plementations of such schemes, since entanglement is eas- 
ily degraded. More recently, a different measure of the 
degree of quantumness of correlations that a bipartite 
system contains has been proposed - the quantum dis- 
cord, which allows an appropriate characterization of the 
quantumness of mixed states [1-7]. This quantity is de- 
fined as the difference between two classically equivalent 
expressions for the mutual information, when applied to 
a quantum system. It has since been shown that systems 
that are not entangled, but exhibit nonzero quantum dis- 
cord, still allow for quantum enhancements in a variety of 
situations [8l4lQj. This suggests that quantum discord is, 
in some scenarios, a more faithful measure of nonclassi- 
cality, as it is sensitive to correlations that are ignored by 
measures like the inseparabilty or the logarithmic nega- 
tivity [5|, [6j. This may allow more robust systems for 
quantum information processing to be recognized. 

Despite the rapidly growing body of theoretical works 
concerning quantum discord, few experiments have in- 
vestigated the topic pj]-[l6|. This may be partly due to 
the fact that analytic procedures to determine discord 
in experimental systems involving mixed states are not 
yet as established as for pure states. Here we measure 
the amount of quantum discord, as well as quantum mu- 
tual information and classical mutual information, that 
are present in twin beams created by the vacuum-seeded 
four- wave mixing (4WM) process in hot Rb vapor. This 
is based on experimental reconstruction of the covariance 
matrix of the system. Additionally, we demonstrate the 
relative robustness of discord as a nonclassicality mea- 
sure, as evidenced by the fact that even for vanishing 
squeezing, non-zero quantum discord still exists. Finally, 
we show that the multi-spatial-mode nature of 4WM re- 



sults in dramatically different behavior of the discord de- 
pending on which pairs of spatial sub-channels of the 
system are used. 

We begin by summarizing the concept of mutual infor- 
mation as applied to both classical and quantum systems, 
and introduce the concept of quantum discord. Classical 
information theory results in two equivalent expressions 
for the mutual information present in a bipartite system 
with two random variables A and £?, which take on the 
values a and 6, with the probabilities pa(o) and £>s(6), 
respectively. The joint probability PAB(a,b) of the sys- 
tem characterizes the correlations between the systems 
A, B. The two alternative formulations for the mutual 
information are |l|: 



I M (A : B) = H(A) + H(B) - H(A, B) 
J(A : B) = H(B) — H(B\A). 



and 



(i) 



Here H(X) = — ^p(x) log 2 p(x) is the marginal en- 
tropy of the probability distribution for a given result X, 
H(A, B) is the joint marginal entropy, and H(B\A) is the 
conditional entropy, describing the average entropies of 
B conditioned on the alternative outcomes of A follow- 
ing B ayes' s rule. Equation (1) determines the amount of 
information one can gain about subsystem A upon mea- 
surement of subsystem B. Quantum discord results from 
the fact that that the two classical formulations of mutual 
information generally do not coincide for a state repre- 
sented by the density matrix pAB- In the quantum case, 
the marginal entropy is replaced by the von Neumann 
entropy S(p) = — Tr(plog 2 p) and yields 



Im(pab) = S(pa) + S(p B ) ~ S(pab)- 



(2) 



The expression for J (A : B) contains the conditional en- 
tropy H(B\A), which, when applied to the quantum case 
results in a term that requires minimization over all pos- 
sible measurements on subsystem A (due to the ambigu- 
ity in defining the state of subsystem B before measuring 
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subsystem A). In general, for the quantum case Im ^ J, 
as here measurements on only one subsystem do not give 
full knowledge about all of the existing correlations. The 
difference is the discord V = Im — J •> given by 

V(A; B) = H(p A ) - H(p B ) + inf ^>^# (p B | fe ), (3) 

k 

where the infimum describes the optimization over all 
possible measurements performed on system A. This 
term makes it an open problem to find an applicable ex- 
pression under many experimental situations. 

To calculate the quantum discord present in a given 
system assumptions about the state of the system must 
be made. We now concentrate on the special case of 
discord for two-mode Gaussian states, for which an an- 
alytic solution for the general case has been given in [E| 
and, specifically for symmetric squeezed thermal states, 
in [l7j. The formalism requires knowledge of the covari- 
ance matrix of the bipartite state. In the following we 
show explicitly how to derive the quantum discord from 
experimentally obtainable data, based on the theoretical 
description of Datta et al. 

A two-mode Gaussian state pab can be described by 
its covariance matrix jab, which may be transformed 
into the standard form [l8|, [l9j with diagonal sub-blocks 
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FIG. 1: Experimental setup for determining the elements of 
the covariance matrix. Two-mode squeezed vacuum is gen- 
erated by a 4WM process. Probe and conjugate beams are 
generated in a cone around the pump beam with an angle 
« 1°. Each beam is separately overlapped with a strong local 
oscillator beam on a 50/50 beam splitter, where both outputs 
are fed into a photodetector. The relative phase of the local 
oscillators is scanned by a piezo mirror in one path. The sum 
and difference of the generated photocurrent signals are fed 
into electronic spectrum analyzers. 5i and 62 are the one- and 
two-photon detunings, respectively. 
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where A = I\ + 

The values n, m, c\ and c 2 are experimentally accessi- 
ble via homodyne and joint-homodyne measurements, as 
discussed in the following. This allows one to express the 
quantum mutual information, the classical correlations 
and the Gaussian quantum discord as a function of the 
symplectic invariants [I?} : 



(6) 
(7) 



Im = KVh) + h(Vh) - h{d+) - h(d-), 



J = h(Vh) - h( x /E™,™), 



VaMiab) = h(y/h) - h(d-) - h(d+) + h(JE™S) 



(8) 



with 



771mm j 

h A\B < 



2ll-{I 1 -AU)(l 2 -\)+2\h\^q-{I 1 -AU)(l 2 -\) 
4(W) 2 



if 



(/i+4/ 4 )(/ 2 + i)/| - 



< 1 



otherwise, 
(9) 



where h(x) = (x + \) log 2 (x + \) — (x — \) log 2 (x — \). 
V > 1 is a sufficient condition for demonstrating the 
entanglement of Gaussian states and V = is reached for 
classical states. A condition arising from the uncertainty 
relation requires d- > 1/2 in order for a state to be 
physical [5[- 

The formulation of discord in Eqs. (8,9) allows for 
direct experimental access via homodyne and dual- 
homodyne detection j20l - [23| . Two- mode squeezed vac- 
uum is generated in multiple spatial modes by a 4WM 
process near the 85 Rb Dl line \5Su>>.,F = 2) |5Pi/ 2 ), 
at A ~ 795 nm as described in [24] and indicated in 
Fig. 1.; probe and conjugate modes are generated at 
±3 GHz relative to the pump. The one-photon detuning 
is Si ^800 MHz and the two-photon detuning 82 is effec- 
tively MHz, as the detected field is determined by the 
local oscillator detuning. Bright local oscillator beams 
are created by seeding a second 4WM process that has 
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FIG. 2: Frequency spectrum for the extracted discord (red cir- 
cles), quantum mutual information Im(pab) (green squares), 
classical correlations J(pab) (black diamonds) and squeezing 
(blue triangles) , respectively (the shot noise limit is indicated 
with a dashed line). These measurements, taken at different 
RF center-frequencies shown, serve to vary the off-diagonal 
covariance matrix elements (that depend on the amount of 
squeezing) , while effectively keeping the on-diagonal elements 
constant (these depend only on the excess noise in each single 
mode) . 



otherwise the same experimental parameters as the first 
4WM process, as in |24j . 

The case examined here involves a symmetric two- 
mode squeezed state which allows for a simple treatment, 
resulting in direct experimental access to the nonzero ma- 
trix elements in Eq. (j4]). A single homodyne detection of 
each of the modes A and B gives the excess noise fluctu- 
ations of the individual modes 

A 2 X A and A 2 X B (10) 

which, when scaled by the shot-noise limit (SNL), re- 
sults in the matrix elements n and m, respectively. The 
off-diagonal matrix elements c\ and C2 are extracted 
from measurements of the fluctuations of the joint- 
amplitude and phase quadratures (difference and sum, 
respectively): 

X_ = (X A - X B )/y/2, Y + = (Y A + Y B )/V2. (11) 

Here Xa, X b and Y A , Y B are the amplitude and phase 
quadratures of states A and £?, respectively. Experimen- 
tally, this corresponds to two separate measurements of 
the sum and difference signals of a joint-homodyne de- 
tection of both modes, while scanning the path length of 
one of the local oscillators relative to the other, as indi- 
cated in Fig. 1. These photocurrents give direct access 
to the variances of the two joint quadratures A 2 X_ and 
A 2 F + when the local oscillator is scanned and the trace 
reaches a minimum. We can then directly determine the 
off-diagonal entries of the covariance matrix 

ci = A 2 X A -A 2 X_ and c 2 = A 2 Y B - A 2 F+, (12) 



FIG. 3: Quantum discord (red circles), inseparability I 
(purple stars), quantum mutual information Im(pab) (green 
squares), classical correlations J(pab) (black diamonds) and 
squeezing (blue triangles, plotted on a linear scale) versus the 
(symmetric) transmission of each two-mode squeezed vacuum 
path when attenuated with neutral density filters. / at each 
point is < 2, confirming that the state is inseparable. The 
theory curves are calculated by using the measured initial 
squeezing and excess noise at a transmission of one. 

after normalizing each quantity to the SNL. 

We experimentally examine the frequency spectrum of 
the discord in our system by changing the RF center- 
frequency at which the measurements are made. The 
joint quadrature squeezing and the excess noise on the 
individual probe and conjugate beams are expected to 
exhibit a somewhat different frequency dependence due 
to the variation of the extra noise from the 4WM process. 
The measured spectra of discord, J&r, J, and squeezing 
are shown in Fig. 2. The noise variance (blue triangles), 
normalized such that 1 corresponds to the SNL and is 
perfect squeezing, exhibits excess noise at low frequencies 
due to technical noise. Optimum squeezing is obtained in 
the range ~ 500 kHz to 1 MHz, with decreasing measured 
squeezing at higher detection frequencies. (All uncertain- 
ties shown are one standard deviation, combined statis- 
tical and systematic uncertainties. The statistical un- 
certainty is obtained from the variance on the measured 
shot noise and excess noise on the single beam measure- 
ments and the joint quadrature measurements. In Fig. 2 
the uncertainties are shown at only one frequency, but 
are representative of typical uncertainties.) Though not 
shown in the plot, there exists nonzero discord for some 
distance after squeezing has reached zero. The frequency 
spectrum of the squeezing is largely reflected in the spec- 
tra of the discord (red dots), the quantum mutual in- 
formation (green squares) and the classical correlations 
(black diamonds), showing that each of these measures 
are equally good for investigating the qualitative varia- 
tion of bipartite correlations. 

The quantum discord captures information about all 
of the elements of the covariance matrix. Another mea- 
sure, the inseparability criterion, defined as / = (AX 2 ) + 
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Transmission 

FIG. 4: Discord as a function of transmission of the local os- 
cillator beams. The local oscillator beams were clipped sym- 
metrically (red crosses) and antisymmetrically (green stars) 
as sketched on the right. The points for « 50% clipping are 
indicated, and the insets have the subregions for 50% clipping 
labelled with Roman numerals. 



(AYj?), only contains information about the off-diagonal 
elements. We now examine how the discord, Jm, J, and 
I scale with the symmetric attenuation of the two-mode 
squeezed vacuum state. Before the double- homodyne 
measurement, we symmetrically vary the amount of at- 
tenuation in both modes of the two-mode squeezed vac- 
uum state with a neutral density filter. This serves to 
change all of the elements of the derived covariance ma- 
trix in a controllable, symmetric manner. In Fig. 3 we 
plot the quantum discord (red dots), the inseparability 
(purple stars), Im and J versus the transmission of the 
state, where we attenuate both beams identically with 
neutral density filters. The state is inseparable for all val- 
ues of transmission measured, according to the criterion 
of / < 2. The discord, J&f , and J are all monotonically 
increasing functions of transmission, and diverge from 
one another as the transmission is increased. The theo- 
retical curves in Fig. 3 are calculated for varying attenu- 
ation (with the standard beam splitter model for losses 
[25l | ) at a transmission of one, with no free parameters. 
While the inseparability and the squeezing (blue trian- 
gles) show linear behavior versus attenuation, the discord 
curve is slightly nonlinear, indicating the influence of the 
on-diagonal elements of Eq. (4). 

The present 4WM system can exhibit multi-spatial- 
mode twin-beam correlations between the probe and con- 
jugate beams [26]. We now investigate the discord and 
quantum mutual information that is present between var- 
ious bimodal sub-channels of the two-mode squeezed vac- 
uum. Our aim is to investigate in this context the sub- 
additivity of information measures and channel capacity 
[271 , [28| . Subadditivity of a channel is denned here as the 
situation when the discord of the total bipartite channel 
is smaller than the sum of the discord between two bipar- 
tite sub-channels. Superadditivity is then the converse. 



Let us define four channels, made up by combinations of 
one of the regions I and II as the left and right halves of 
the probe mode, and one of the regions III and IV as the 
left and right halves of the conjugate mode (as shown in 
Fig. 4). The correlations between spatial regions in the 
probe mode and the conjugate mode are expected to be 
centro-symmetric with respect to the pump beam, as in- 
dicated in Fig. 1, which allows us to differentiate between 
mostly correlated and largely uncorrelated subregions. 

We measure the quantum discord, V{II\IV), present 
between modes II and IV by cutting the local oscillator 
beams with a knife-edge, and only allowing these modes 
to be fed into the dual-homodyne detection system. This 
serves to project out only modes II and IV. As seen in 
Fig. 4, the quantum discord present between these two 
channels, V(II; IV), is much less than half of the total 
discord present between the entire, uncut probe and con- 
jugate modes, V(I + II; III + IV). In this case, the 
discord is super additive, or in other words, we have cho- 
sen subregions of the beams that are largely uncorrelated. 
(The same result is obtained with V(I;III), confirming 
the spatial symmetry when clipping the beams in this 
direction). In contrast, the discord between modes II 
and III, V{II] III), is significantly more than half of the 
total discord between the total, uncut modes. In this 
scenario, the discord behaves subadditively. This subad- 
ditive property (2?J of the discord may allow for more 
effectively exploiting the channel capacity of a total sys- 
tem if it is correctly divided into various bimodal sub- 
channels. This shows that the discord varies dramat- 
ically in multi-spatial- mode squeezed vacuum systems, 
depending on which spatial channels are utilized. Here, 
when symmetric channels are used, the discord is subad- 
ditive with respect to the unaltered system. We also find 
that when asymmetric channels are chosen, the amount 
of discord quickly decreases, resulting in super additive 
behavior. 

Existing theoretical work concerning quantum discord 
indicates that discord can in principal be employed in 
various quantum computation and quantum communica- 
tion schemes to provide enhancements beyond those al- 
lowed classically. In particular, many schemes require op- 
tical delay and storage of quantum correlations that are 
always accompanied by some loss and decoherence. The 
relative robustness of quantum discord with respect to 
decoherence effects provides motivation to more closely 
investigate these systems experimentally. Specifically, 
the present scheme for the determination of discord is 
applicable if one is interested in the evolution of quan- 
tum correlations in systems that involve stored light [29| 
and fast light [30|, HH , where squeezing and entanglement 
may be easily degraded. Discord as a measure of quan- 
tum correlations in mixed bipartite states may offer a 
means to experimentally test the non-broadcasting theo- 
rem for the class of quantum-correlated bipartite states as 
a simple example of mixed states [32|, [33| under fast-light 
conditions, and how classical and discordant states can 
be distinguished in this regime. Examining the amount of 
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discord, rather than squeezing or entanglement, present 
under these exotic conditions should allow for a more 
easily experimentally obtainable measure of any "quan- 
tumness" that is present in the system. 
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